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TURBULENT BOUNDARY -LAIER TEMPERATURE RECOVER! FACTORS 
IN TWO-DIMENSIONAL SUPERSONIC FLOW 
By Maurice Tucker and Stephen H. Maslen . 

SUMMARY 

An analytical method is presented for obtaining the turbulent 
temperature recovery factors for a thermally insulated surface in 
supersonic flow. The method is an extension of SqiAire's analysis 
for incompressible flow. The boundary -layer velocity profile is 
represented by the power law and a similarity is postulated for the 
squared -velocity and static -temperature-difference profiles. 

The emalysis indicated that the recovery factor decreased with 
increasing Mach number. For the range of Prandtl numbers considered 
(0.65 to 0.75), the recovery factors at a stream Mach number of 10 
were, on the average, about 5 percent lower than the limiting values 
at zero Mach number. The Reynolds number effect on recovery factor 
was of secondary importance. An approximation formula that represents 
the computations to within 1 percent is Included for engineering 
calculations. 


INTRODUCTION 

Detennination of the temperature attained at the surface of a 
thermally insulated plate in lamlmr boundary -layer flow in the 
absence of radiation (a particular case of the plate -thermometer 
problem) has been the subject of several analytic investigations. 

The surface temperature Tav ms-y be presented in terms of the stream 
static temperature tq, the temperature recovery factor r, and the 
local stream Mach number parameter m^ as 


■^aw “ ^l(l + rm^) 


( 1 ) 


(All symbols used in this report are defined in appendix A.) For . 
very low speed flows, Pohlhausen (reference 1, pp. 627 -631) found 
that r vas a function only of the laminar Prandtl number Pr and 

could be approximated as r = Pr^/^. References 2 and 3 indicate 
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that the Pohlhausen approximation is also .applicable to supersonic 
laminar boundary- layer flow. 

Existing analyses of the plate -thermometer problem for turbu- 
lent boundary-layer flow (see, for example, reference 4) are limited 
to the case of constant fliiid properties* The results of refer- 
ences 5 and 6 indicate that, for the low-speed flows considered, the 
turbulent recovery factor may be represented by the approximation 
j. _ pr^' The analysis conducted at the NACA Lewis laboratory and 
presented herein essentially extends the method of reference 6 to 
supersonic two-dimensional turbulent boundary-layer flow* 


ANALYSIS 


Development 

The principle of conservation of energy req.uires the total- 
energy increment transported throu^ any section normal to a thermally 
insvilated surface along which a two-dimensional boundary -layer flow 
exists to be constant; that is. 



pucp(T-Ti) dy 


= constant 


( 2 ) 


^ere the specific h^t Cp is taken as a constant. Inasmuch as no 
energy transport occurs at the leading edge of the plate, the con- 
stant appearing in equation (1) must be identically zero. Squire'- s 
basic equation (reference 6) for obtaining the temperature recovery 
factor is then given in the present notation as 


Pu 

Pl"l 


c (T-Ti) dy = 0 


0 


(3) 


The treatment of reference 6 is confined to the case of low-speed 
flow by assuming that P = Pi* The method presented herein does not 

contain this limitation* 

For fluids with Prandtl number less than 1, the thermal boundary- 
layer thickness A exceeds the dynamic boundary-layer thickness 8. 
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For a Prandtl number of 1, in which case A » 6, the following static- 
temperature profile (from equation (7) of reference 7) is obtained; 


- ti 


= 1 + 


s 1 + m^ j^l - g^(y)J 


(4) 


which Indicates a similarity of the squEired -velocity and static- 
temperature -difference profiles. Squire (reference 6) makes the 
assumption that such similarity also exists for fluids with Prandtl 
number differing from 1. The assumption appears plausible when the 
Prandtl number is relatively constant throu^ the boundary layer 
and does not differ too greatly from 1. 


In the analysis^ the turbulent Prandtl nvunber 


is 


(w)'T> 

assumed to be constant along any boundary -layer section and equal 
in value to the laminar Prandtl number iaCp/k at the surface. These 
assunrptions are consistent with the conclusions of references 8 to 10, 
namely, tha,t the turbvilent Prandtl nxauber has a constant value of 
about 0.7 regardless of the la,mjnar value. Although these assumptions 
restrict the analysis to a consideration of fluids having a laminar 
Prandtl nmber near 0.7, the restriction is not a serious one inasmuch 
as the Prandtl nvunbers of most gases are in this range. 

The slatic-temperattire-difference profile for A 8 is assumed 
to be related to the profile for A « 8 by a constant scaling factor 
Tj = a/ 8 and the recovery factor r. For compatability with equa- 
tions (1) and (4), the static -tenperature profile for fluids with 
Prandtl number different from 1 is then written as 


t 1 -,~2 
— = 1 + rm 
tl 


(0 


rui 


1 + 


2cpti 


1-g' 


(0 


(5) 


This temperature relation, postulated by Squire (reference 6), will 
'be used in the analysis. The following expresssion is thereby obtained 
for the quantity Cp(T-Tj^) ; 


u2r 

)(T-Tl) = 


1-g 


(O] 


( 6 ) 
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When the perfect gas law euid equations (5) and (6) are oised^ equa- 
tion (3) takes the form 


0 





^ g(y) dy 

— Z-f 

1 rm 

‘-'(f)] 



(7) 


The assumption is now made that the velocity profile can he repre- 
sented hy the power law 


i/w 


( 8 ) 


The quantity q 


s is defined and 

g(y) = P 

for 

0 < y< 6 

g(y) = 1 

for 

y ^ 8 


for 

0 < y s A 

4 ^) - " 

for 

y. ^ A 

can he written as' 


r 


. -1 pi 

T]^/^ (r-l) + 

fr,2/N 

1 

> 

% 

V 

' JJo 

r 


Tir / 


m 


^vl/N A-q2 


(0 


(9) 


where A E 


_ 1 + rm^ 


rm 
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As is shovn in appendix B# the scaling factor ti Is related to 
toe recovery factor r "by 

r « n^Pr^ » 

This relation \«as also obtained for incompressible flow in refer- 
ence 6. The effect of a variable scaling factor on toe recovery- 
factor was, also considered. By .use of the conditions of reference 6 
(p = and A < 5)/ computations were made for toe scaling fac- 
•bor S given by 

(l-D/l / Wi< 

s » n pt] < I < 


s « p 

For variable scaling factors, equation (lO) takes the form 

r = S^(y^)Pr aoa) 

When it is assumed -that L » r and P = 0.2, toe recovery factor 
ob-tained is 98.6 percent of toe value obtained for constant scaling 
factor S * T]. In view of toe relative insensitivity of toe computed 
recovery factor to relatively large changes in scaling factor for 
■values of N -that sire appropriate to turbulent boundary -layer velocity 
profiles, use of. a constant scaling factor was considered adeqmte. 

When toe substltutlcne 

and 

a1/2 

q » A ' z 


1/L* y 

n 
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are ina.de and eq^uation (10) Is used^ equation (9) takes the form 


^AB 


i/n iW 




■) (aI/V/®") 


B- 


0 


z. dz 

l-z2 


-A- 


m 


- 1/2 - 1/2K 

A B 


z*'-^ dz 
l-z2 


^(2-N)/2 ^(2-N)/2N 
N+I 


(2-N)/2 - 1/2 (2-N)/2 

rA B xk 

N(N+l) • N 


F = 0 


( 11 ) 


The ten^ierature recoyery factor r as a function of Mach’ number para- 
meter m, velocity -profile parameter N, and laminar Prandtl number Pr 
is ob-tained from equation (ll). The effect of chord’wlse pressure 
gradient on recovery factor enters only to the extent 'that "the velo- 
city profile is dependent upon -the gradient. 


Me-thods of Solution 

An explicit solution of equation (ll) for -the recovery factor r 
is obviously impossible inasmuch as r appears in both the quan- 
tities A and B. The following procedure -was employed: A me-thod •was 

de-vlsed for obtaining a first approximation to "the recovery factor. 
Values ob’tained by ■this me’thod were -then improved by use of the Ne'wton- 
Baphson iteration me’thod (reference 11). De’tails of the procedure 
follow: 

* . 

First approximation me’thod. - For Prandtl numbers near 1, it is 
expedient ’to ’write the recovery factor as 

r = Pr®" =[^1-(1-Pr)] s (1-H)®’ 

and to approximate the quantity (l-H)°' as 1-oH. The integrals 
appearing in equation (ll) may be denoted by I « I(m,N,H). This 


( 12 ) 
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functional relation 
witii the assumption 

l(m,N,H) = l(m,W,0) 
these substitutions 


may he expanded in a Taylor's series, 
that the Prandtl number is near 1, as 


consistent 


5hJ 


. " H *=0 

in eqmtion (ll) yields 


Introducing 


ai 


1 - (N+1)Q 


1 - N + (N+1)Q 


N -1-m^ 
1+m^ 


where 


Q 




(H+l)/2 




- 1/2 


= L 


dz 

l-z^ 




(13) 


For application of equation (12) , punch-card equipment ms used to 


evaluate the definite integral 
from 4 through 11. 




*2 for integer values of 


1-z 


|0 



Iteration method . - With an approximate solution ra of equa- 
tion (ll) obtained from the method just described, an improved value 
is given by the Newton-Raphson method as 


F(ra) 

^b = ^a - 

• . 


The quantity F is defined by equation (ll) and the partial deri- 
vative is obtained as 
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A- 

2Nr 



- 1/2 - 1/2N 
PA B 



dz 

l-z2 


+ 


A- 

2N(w+1) 


2{N+1) + B^^^^A+N(l-A)-2j - | ^A+N(1-A)J B 


(2-N)/2N 


(15) 


Procedure . - First approximations for r using equations (12) 
and (13) were obtained for Prandtl nuDibers of 0.65, 0.70, eind 0.75 
at local stream Mach numbers of 3, 6, and 10 with •values of -the pro- 
file parameter N of 5, 7,- 9, eind 11. Two iterations, oh'tained hy 
use of equation (14) , were vised for each case to give results correct 
•to four decimal places. 


Limiting solutions. - The limiting solutions of equation ( 11) 
are of interest. For "the limiting case of infinite Mach number, the 
recovery factor rj^ is ob'tained as 


r 

CD 


B 


,1/a 




- 1/2N 


z^ dz 

l-z2 


_ 1 

B 2 1 

N(N+1) " N 


g(2-N)/2N 

N+1 


(16) 


Compu'tations indicate -that r is practically independent of "the 

Oo . 

profile parameter N and can be represented wi'fchin 1 percent by 
the relation 

r^ » 0.670 Pr + 0.322 (17) 

for Prandtl numbers in -the range from 0.65 to 1. 

For -the limiting case of zero Mach number, expansion of ‘the 
integrands in equation (11) permits writing "the equation as 
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N(N+1) 


- 1/2 3/2 

2Pr 

W(N+2) 


i/n (n-i)/m 
. '"0 2 

CN+2)(H+3) CN+l)(B+3) 


2 

Approximating Pr*' » (1-H)^ as 1-GH + and solving for rQ 

by the Kewton-Eaphson metbiod gives the result 


Pr 


N + 1 
3N + 1 


(18) 


The expression for rQ is Identical vlth Squire's result (refer- 
ence 6) for incompressible flov of a fluid with Prandtl number 
greater than 1. The limiting values of recovery factor obtained by 

N + 1 

the first approximation method are rQ = Pr^^ ^ ^ and = Pr. 

The approximation method thus overestimates the effect of Mach 
number on recovery factor. 


RESULTS AMD DISCUSSION 

The coniputed recovery factors are given to three decimal places 
in table T. It may be noted here that the first approximations 
obtained by use of equations (12) and (13) were found to be within 
1.4 percent of the values tabulated. The maximum error occurred at 
a local stream Mach number of 10 for a Prandtl number of 0.65 and a 
profile parameter N of 5. The following expression 

N + 1 + 0.528 
3N + 1 + 

r - Pr 

represents the tabulations to within 1 percent and may. prove of con- 
venience in engineering calculations. With the recovery factor 
known, the stagnation temperature profile for constant scaling factor 
may be obtained from equation (6) as ' 
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A typical profile is shown in figure 1. The discontinuities in slope 
at y/S = 1 and y/5 = t] eind the incorrect value of the slope at 
the wall result from xise of a pover-law profile throio^out the houn- 
dary layer. 


Effect of ^ch Number 

Table I indicates that the turbulent temperature recovery factor 
decreases with an increase in local stream Mach number. For the range 
of Prandtl numbers and velocity profile parameters considered, the 
decrease is about 5 percent for a change in Mach number from 0 to 10. 

A typical variation of turbulent recovery factor with Mach number is 
shown in figure 2. The corresponding laminar recovery factor obtained 
by the method of reference 2 is also shown. The invariance of laminar 
recovery factor with change in Mach number apparently results from use 
of the assumption u ~t. The results of reference 12 indicate that 
use of the assumption ja ~ t*^*® leads to a decrease of lewlnm- recovery 
factor with increase in stream ^fech number. For a Prandtl number of 
0.70, the recovery factor at a nominal Mach number of 5.4 is 0.815 as 
compared with the value 0.837 for incompressible flow. 


Effect of Reynolds number 

It is well known (reference 1, p. 340) .that use of the power law 
(equation (8)) to represent the turbulent boxaidary -layer velocity pro- 
file requires that the parameter N increase with Reynolds number. 
Equation (4b) of reference 7 presents the following approximate guide 
for such variation; 


l/l4 


s,w 


N « 2.6 R 


( 21 ) 
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The emalysis herein, as do low-speed flow analyses, thus req^uires the 
recovery factor to increase with increasing Reynolds number for a 
given Mach number. This increase is greater at high Mach numbers. 

In general, however, the effect of Reynolds nmber upon recovery 
factor is, according to the analysis, a secondary one. From the 
relative insensitivity of the results to changes in profile para- 
meter N, the effect of moderate chordwise pressure gradients upon 
the recovery factor should also prove of secondary importance. 


Comparison With Experiment 

Experimental turbulent recovery factors obtained at a ^fe.ch 
number of 2.4 are presented in reference 13 for steady flow over a 
flat plate with natural transition in a tunnel having test-section 
1 1 

dimensions of 52 by 5^ inches. As the Reynolds number based on dis- 
tance from the plate leading' edge and on free-stream kinematic vis- 
cosity was increased from 2 x 106 to 6.7 x 10®, the turbulent recovery 
factors decreased from 0.897 to 0.884. The trend with Reynolds 
number is thus in the opposite direction to that predicted by the 
various analyses. 

In reference 13, the Prandtl number at the plate surface was 
taken as 0.72 and the profile parameter N as 5. Using equation (19) 
as an interpolation form'ola and substituting these values of Prandtl 
nvtmber and profile parameter yields, at a Mach nmber of 2.4, a 
recovery factor of 0.872, which is from 1.4 to 2.8 percent lower 
than the extremes of the experimental -values. The profile -pairameter 
•value N = 5 chosen in reference 13 appears somewhat low for the 
range of Reynolds numbers used. A -value of N = 7 ob-fcained from equa- 
tion (2l) by averaging -with respect to plate -stirface and s-tream Rey- 
nolds number yields a recovery factor of 0.877 for Pr = 0.72 and 
M = 2.4. This value of recovery factor is from 0.8 to 2.2 percent 
lower -than -the extreme values measured. The -theoretical value of 
laminar recovery factor given in reference 2 for Pr = 0.72 is 0.845 
as compared -with -the experiraen-tal -value of 0.881 given in reference 13. 
More extensive tests, particularly at very high Mach numbers, are 
required for a decisive check of -the assumptions involved in the 
analysis. 
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SIM4AEY OF EESULTS 

The present work essentially extends the low-speed analysis of 
Squire to the case of two-dimensional supersonic flow along a ther- 
mally insulated surface. The analysis makes use of the following 
simplified model of turbulent boundary -layer flow: 

The turbulent Prandtl ntunber was assumed to be constant along 
any bomdary -layer, section and equal in value to the laminar Prandtl 
number at the svirface. The boundary -layer velocity profile was 
approximated by a power law. The presence of the laminar sublayer 
was not considered except for evaluation of the surface Pra n dtl 
number and the similarity of the squared -velocity and the static- 
temperature -difference profiles postulated by Squire was assumed. 

Under these conditions, the temperature recovery factor in tur- 
bulent flow was shown to decrease with increasing Mach number. For 
the range of Prandtl numbers considered (0.65 to 0.75), the recovery 
factor at a Mach ntmber of 10 was about 5 percent lower than the 
limiting value at a Mach number of 0. As in low-speed analyses, the 
recovery factor was shown to increase with an incr^se in Reynolds 
nxmiber. The Reynolds number effect upon recovery factor was, in 
general, of secondary Importance. The approximation for recovery 

N + 1 + 0.528 
3N + 1 + 

factor r = Pr (where Pr is the laminar Prandtl 

nvimber, N is the profile parameter, and is the local stream 

Mach number) was found to represent the computed results to within 
1 percent. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, October 27, 1950 
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APPEMDU A 
Synibols 

The following synibols are used in this report: 

l+na^ 

rm 


B S 


F 

a 

s(y) 

H S 
k 

L • 

M 

n .2 

m 

N 

Pr 

Prt 

P - 
<1 = 


_r_ 

Pr 


specific heat at constant pressure 
specific heat at constant volume 
defined by eqmtion (ll) 
generalized exponent 
velocity-ratio function, u/u^ 

1 - Pr 

thermal conductivity for laminar flow 
limit of integration 
Mach number 

Mach number parameter, 21^ M^^ 

velocity- profile parameter, u/uj^ - 

hCp 

Prandtl number for Inmlnnr flow, — ^ 
Prandtl number for turbtilent flow, 

iCt p 

i/n 
i/n 


i/n 


(i) 
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E 


s,w 


Reynolds nvimber "based on distance from leading edge of plate 
flnri stream or wall value of kinematic viscosity 


T -t, 
aw 1 


temperature recovery factor, 


S 

T 

t 


u,v 


x,y 


z 

a 

3 

7 

A 

6 

c 

n 

p 


scaling factor 
stagna.tion temperature 
static temperature 

velocity components parallel and normal to surface, respectively 

cartesian coordinates parallel and normal to surface, 
respectively 

variable of integration 

recovery-factor exponent, R = Pr°’ 

eddy conductivity 

ratio of specific heats, Cp/cy 

thermal houndary -layer thickness 

dynamic, houndary -layer thickness 

eddy viscosity 

scaling factor, a/6 

laminar coefficient of viscosity 

density 


Subscripts : 

1 local stream 

aw at surface of adiabatic or thermally insulated plate 
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0 at limiting case 

oc at limiting case 


Ml = 0 
Ml = “ 
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APPENDIX B 


Eelation Between Scaling Factor i] and Becovery Factor r 

The following equation for the conservation of energy Is assumed 
to hold for steady two-dimensional compressible turbulent boundary- 
layer flow; 


dt Bt ^ / X 

PUOp gj + PTOp 3^ - U S 


f 3u\ ^ ^ F/v 


3t 

¥ 


(Bl) 


Adding to this equation the corresponding equation of motion multi- 
plied throu^ by the velocity u yields, after some manipulation. 


Pu 


6t b ( v?\ 


6t 6 /u 

A] 

°p 55E ^ \T) 

+ Pv 

°P 5? ^ 





(p«) 


3? * 5? (rjj ^ ^ I (pK)op ^ 


(B2) 


In deriving equation (B2), the assumption has been made that the t\xr- 
bulent Brand tl number Prt = cp is constant. At the plate 

surface where u = v = 0, equation (B2) can be written as 




pT'-ipTl) 






or, noting that 



=0, then 

y=0 


(h+e) 


(PpT-PpTl) 


y=o 


(l-Pr^) (h-k) 

y*0 



(B3) 
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The following relation may be obtained from eqmtion (6) ; 


[ 


^y2 


y=0 


H) f { 


du 


y*0 


(B4) 


The desired relation between n and r is obtained from equa- 
tions (B3) and (B4) as 


r = Pr^ = Pr 


( 10 ) 


in view of the assumption regarding the equality of the turbiilent 
and laminar Prandtl numbers at the surface (y=0)» 
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TABLE I - VABIATION OF TURBULENT TEMPERATURE RECOVERY 
FACTOR WITH PRAKDTL NUMBER Pr, VELOCETY- 
PROFILB PARAMETER N^ AND STREAM 
MACH NUMBER M^ 


Prandtl 

ntmiber 

Pr 

Velocity- 
profile 
parameter, N 

Mach number, M^^ 

0 

3 

6 

10 

0.65 

5 

0.851 

0.834 

0.813 

0.796 


7 

.855 

.840 

.821 

.803 


9 

.857 

.844 


.808 


11 

.859 

.847 

.830 

.813 

0.70 

'5 

0.875 

0.860 

0.842 

0.826 


7 

.878 

.865 

.848 

.833 


9 

.880 

.869 

.853 

.838 


11 

.882 

.872 

.856 

;841 

0.75 

5 

0.898 

0.885 

0.870 

0.856 


7 

.901 

.890 

.875 

.862 


9 

.902 

.893 

.879 

.866 


11 

.904 

.895 

.882 

.869 




Normal distance parameter, y/fi 


20 


MCA TN 2296 



Figure 1, - Stagnation temperature profile for turbulent flow at a stream Mach number of 3 for 
Prandtl number of 0*70 and boundary-layer velocity profile parameter N = 7,. 
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Figure 2. - Variation of temperature recovery factor with stream 
Mach number In two-dimensional flow of fluid with Prandtl number 
of 0.7. For turbulent flow, boundary -layer profile parameter 


N = 7. 
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